Tunnel Magneto-resistance in GaMnAs: going beyond Julliere formula. 
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The relation between tunnel magneto-resistance (TMR) and spin polarization is explored for 
GaMnAs/GaAlAs/GaMnAs structures where the carriers experience strong spin-orbit interactions. 
TMR is calculated using Landauer approach. The materials are described in the 6 band k ■ p 
model which includes spin orbit interaction. Ferromagnetism is described in the virtual crystal 
mean field approximations. Our results indicate that TMR is a function of of spin polarization and 
barrier thickness. As a result of the stong spin orbit interactions, TMR also depends on the the 
angle between current flow direction and the electrode magnetization. These results compromise 
the validity of Julliere formula. 

PACS numbers: 



The relative orientation of the magnetization of two 
ferromagnetic electrodes can affect dramatically electron 
transport across a tunnelling barrier connecting them. 
This phenomenon gives rise to the so called tunncl- 
magneto-resistancei (TMR), TMR = {R AP - R P )/R AP 
where Rap and Rp are the resistances for anti-parallel 
and parallel orientations respectively. TMR is exploited 
to fabricate devices which are ultra-sensitive to variations 
of an external magnetic fieldS. Microscopic understand- 
ing of TMR, based upon the hypothesis that spin is con- 
served in the tunnelling process, leads to the well known 
Julliere Formula 3 - 4 - 



TMR 



2P L Pr 
1 + PlPr 



(1) 



which relates TMR with Pl,r, the polarizations of the 
left (L) and right (R) electrodes. Assuming Eq. JTJ 
is correct, it permits to extract the spin polarization of 
the electrodes from the value of TMR in symmetric tun- 
nel junctions, regardless of the physical properties of the 
barrier. Eq. Q can be derived doing second order per- 
turbation theory in tunnelling amplitude 5 . 

Diluted magnetic semiconductors like Gai-^Mn^As 
are ferromagnetic below Tc — 150 K&. This type of 
materials rise a lot of interest because they afford the 
integration of ferromagnetic and semiconducting func- 
tionalities in a single device^. Substitutional impu- 
rities of Mn in GaAs are acceptors. The holes re- 
leased by the Mn are responsible of the magnetic or- 
dering and transport properties of Gai-^Mn^As. Sev- 
eral groups have been able to fabricate tunnel junctions 



* brey@icmm.csic.es 



with GaMnAs in the electrodes2ii2iii and they have re- 
ported large values of TMR finding that its value de- 
pends on the properties of the barrier 9 . Furthermore, 
spin orbit (SO) interactions for the holes in the valence 
band of GaAs is very strong (Ago — 0-34 eV) so that 
the spin of the holes is not a good quantum number. 
Both the experimental results2*i2*ii and the failure of Jul- 
liere's hypothesis on spin conservation, lead us to study 
TMR in GaMnAs/AlAs/GaMnAs systems using a non- 
pertubative approach which fully includes spin orbit in- 
teractions. We anticipate our main result: in GaMnAs 
based heterostructures TMR depends on both the bar- 
rier thickness, (4, the angle formed by the current flow 
and the magnetization. All these features depend on the 
strength of the SO coupling and are relevant for the use of 
GaMnAs based heterostructures in magnetoelectronics. 

In this letter we calculate vertical trans- 
port (z-direction) in epitaxially grown GaM- 
nAs/GaAlAs/GaMnAs tunnel systems. We take the 
same Mn density at the two electrodes. The materials 
are described in the 6 bands k • p model which captures 
the interplay of SO interaction and the threefold orbital 
degeneracy in the top of the valence bands. Ferromag- 
netism is originated by the exchange interaction between 
the spin of the localized Mn and the itinerant holes. In 
the case of metallic samples the mean field and virtual 
crystal (MF-VC) approximationiiiii 5 ^^ accounts for 
a number of experimental facts like the dependence of 
T c on both the Mn and hole densities^, the magnetic 
anisotropyia and magnetic circular dichroism 1 — . In this 
approach spontaneous magnetization is characterized by 
effective Zeeman magnetic field Ti. which represents the 
effect of exchange interaction of the Mn magnetization 
on the spin of the holes. The MF-VC approximation 
restores the translational invariance in the Hamiltonian, 
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FIG. 1: TMR Z as function of P. d b =15A and Vf,=300meV. 
Case (a)p=0.1nm~ 3 and (b)p=0.4nm~ 3 . 

and it is possible to label the eigenvectors and eigenval- 
ues by a band index ,n, and a wavevector (fc z ,kii). Due 
to the SO interaction, energy bands and Bloch states 
depend on both the magnitude and the orientation of Ti, 
and spin is not conserved. For a given set of parameters 
representing GaMnAs, a magnetic orientation SI and a 
temperature, the above approach yields the carrier spin 
polarization P. In this model, the electronic structure of 
GaMnAs is fully characterized by P, the hole density p 
and il. For a fixed Mn-hole exchange coupling strength, 
p and SI there is a one to one correspondence between 
the Mn concentration and P. We describe the valence 
bands of the barrier (GaAlAs) with the same parameters 
used for GaAs, but offse t 17 i 18 i 19 a potential V& of the 
order of 300meV. 

The conductance is calculated using the Landaucr 
formula 17 ! 18 . 19 , 

G ti L ,n R = XI T u!L'( e f) ( 2 ) 

n,n' ,ky 

where Ep is the Fermi energy, T n ^ n {Ep) is the trans- 
mission probability from a state in the left with band 
index n to a state in the right with band index n' , 
and Ql.r are the magnetization orientation of the elec- 
trodes. Only the energy and k|| are conserved in the 
tunnelling process. The transmission matrix is obtained 
in the transfer matrix method applied to the 6 bands kp 
hamiltonianiZii^ii 9 . including the exchange field 

Conventionally, G depends on the angle formed by Qj, 
and f2fl. In presence of SO coupling the conductance 
depends also on the angle formed by the current with the 



FIG. 2: TMR Z and TMR X as function of d b for P=0.75 and 
P=0.89. V b =300meV and p=0.1nm" 3 . 

magnetization. This makes it necessary to define both 
TMR Z , for the case of current parallel (or antiparallel) 
with the magnetization, and TMR^, for the case where 
the magnetizations are perpendicular to the current flow: 

Gt t — Gt i G^_ <_ — G« > 

TMR Z = ^ TMR X = — ^= . 

(3) 

where we denote the positive (negative) z-direction as 
t (|), and the positive (negative) x-direction as (— >). 

In Fig.l we plot TMR Z for a symmetric GaM- 
nAs/ AlAs/GaMnAs junction as function of the bulk spin 
polarization P for two hole densities p=0.1nm -3 and 
p=0.4nm -3 . The barrier height is Vb =300 meV . Both 
the Julliere expression and TMR Z are increasing func- 
tions of the electrode spin polarization P. However for in- 
termediate values of the polarization and p=0.4nm the 
difference between the Julliere and the calculated TMR 
can be of a factor of two. The discrepancy remains if we 
use the Fermi surface polarization^ P FS instead of the 
bulk polarization P . In order to address whether the 
discrepancy is due to the SO interaction, we also show 
in figure 1 TMR Z in the case where the SO interaction is 
suppressed in the barrier. Interestingly, TMR Z is larger 
when SO is suppressed and closer to the Julliere result, 
but still off. 

In Fig. 2 we plot both TMR Z and TMR X as a func- 
tion of the barrier thickness db for two different values of 
P. The Julliere value, independent of db is also shown. 
Due to the strong SO coupling in the system, TMR 2 
and TMR X can be quite different, and because of the 
complicated matching at the interface, its relative mag- 
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FIG. 3: Anisotropic TMR as function of the bulk polarization. 
The error bar is an estimation of the numerical error in the 
calculations. 



and spin polarization of the electrodes. As in the conven- 
tional case, our results show that TMR is an increasing 
function of the electrodes polarization P, quantitatively 
different than Eq. (JIJ. Large values of TMR can be 
reached for sufficiently high P even in the presence of 
strong spin orbit interaction. In contrast with Eq. ifTjl. 
TMR depends on the barrier thickness. Therefore, JJJ is 
not an appropriated tool for inferring the spin polariza- 
tion of the electrodes in GaMnAs based heterostrucures. 
Finally, because of the spin orbit interaction, both the 
TMR and the conductance depend on the angle between 
the current and the magnetization orientation. These 
are qualitatively new physical phenomena that might be 
exploited to build spin valves with new functionalities. 
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nitudc depends strongly on the parameters of the sys- 
tem. In agreement with experiments^, and tight-binding 
calculations^ we find that TMR decreases rapidly for 
thin barriers, in marked contrast with Eq. Q). We also 
plot the TMR when the magnetization is oriented in the 
[101] direction, TMR ZX . This quantity turns out to be al- 
most exactly the average of TMR X and TMR Z , indicating 
that the experimental difference in the TMR between 
the [100] and [110] field directions is due to the SO cou- 
pling and not to the cubic magnetocrystalline anisotropy 
in the GaAs crystal structure. 

A quantitave comparison between our results and the 
experiments is difficult, because quantities like the carrier 
density and band offset are unknown in the heterostruc- 
tures. Even db is not accurately knownS. Qualitatively, 
though, our results account for the decreases of TMR 
with the barrier thickness^, and also show a strong de- 
pendence of the TMR values on the angle formed by the 
magnetization and the current flow. 

Impurity scattering combined with the strong SO cou- 
pling produces anisotropy in the DC transport proper- 
ties of bulk GaMnAsSi. We have analyzed the variation 
of the ballistic tunnelling resistance when the current 
flow is parallel or perpendicular to the magnetization. 
In Fig. 3 we plot the anisotropy tunnelling magneto resis- 
tance (ATMR), 

ATMR = GtJ ~ G ^'^ . (4) 

The tunnelling current can change up to 6 % for large 
values of P. This number is comparable with those ob- 
tained for bulk 21 . Therefore, in such a heterostructure, 
the ATMR can add a significant contribution to the bulk 
anisotropy magneto resistance of the electrodes. 

In conclusion, the presence of strong spin-orbit inter- 
action modifies the conventional relation between TMR 
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Figure Captions 



Figure 1: TMR Z as function of P. <4=15A and 
H=300meV. Case (a)p=0.1n m - 3 and (b)p=0.4n m - 3 . 



Figure 2: TMR Z and TMRj as function of d b for P=0.75 
and P=0.89. V 6 =300mcV and p=0.1n m - 3 . 



Figure 3: Anisotropic TMR as function of the bulk polar- 
ization. The error bar is an estimation of the numerical 
error in the calculations. 



